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The equations for waves on the surface of an irrotational incompressible fluid are derived in the 
coordinates of the velocity potential/stream function. The low frequency shallow water approxima- 
tion for these waves is derived for a varying bottom topography. Most importantly, the conserved 
norm for the surface waves is derived, important for quantisation of these waves and their use in 
analog models for black holes. 
PACS: 47.90.+a, 92.60.Dj, 04.80.y. 



CN . I. INTRODUCTION 

^ ■ 

o ■ 

ps) ' One of the most fascinating predictions of Einstein's theory of general relativity is the potential existence of black 
holes - i.e. space-time regions from which nothing is able to escape. Perhaps no less interesting are their antonyms: 
' white holes which nothing can penetrate. Both are described by solutions of the Einstein equations and are related 
, to each other via time-inversion, see e.g. [H, 

' It is equally fascinating that some of the predictions for fields in a black hole spacetime can be modelled by waves 
' in a variety of other situations, with the interior of the black hole or white hole horizons can be mimicked by fluid flow 
which exceeds the velocity of the waves in some regions. One of these is the use of surface waves on a incompressible 

fluid [3]. One can alter the flow properties of the fluid by placing obstacles into the bottom of a flume (a long tank 

' along which the water flows) to speed up and slow down the fluid over these obstacles. 

Q^' One of the difficulties in the theoretical treatment of such systems is the complicated boundary conditions on the 
bottom of the tank (where the fluid velocities must be tangential to the bottom) and the top (where the pressure of 
the fluid must be zero or at a constant atmospheric constant pressure). In fact as we will see the equations for the 
fluid itself are remarkably simple. The interesting physics arises entirely from those boundary condition. 

We will be interested in irrotational, incompressible flow. While both are certainly approximations for water flow 
(the former assumes no turbulence, and no viscosity which would create vorticity at the shear layer along the bottom, 
, while the latter assumes that the velocity of sound in the fluid is far higher than any other velocities in the problem) . 
' While this problem has been investigated before this is in general in the three dimensional context (which is more 
C J difficult) and using approximations and expansions for the shape of the bottom. 

I will assume that the fluid flow is a two dimensional flow- ie is uniform across the tank and that the tank maintains 
[ a constant width throughout. This is much simpler case than three dimensional flow, which allows the coordinate 
. transformations I use. 

■ The usual spatial coordinates are x, y with x being the horizontal direction in which the fluid flows, and y is the 
vertical direction (parallel to the gravitational acceleration, g, directed in the negative y direction. 
The Euler-Lagrange equations are 

dtv + V ■ \7v = —gsy — V- (1) 
H ■ P 

V-v^O (2) 

where the second equation is the incompressibility condition. In the usual way, if we assume that the flow is irrotational, 
then 



X 



v = W4> (3) 



And the above equation can be written as 



y {dtcP+^v^+9y + ^^ =0 (4) 
V^4> = (5) 

_ £ 

p 

In the following I will consider only flows in the x — y directions. Everything is assumed to be independent of z. 
Consider the vector w — Cz x v. This vector also obeys 

V • zi? = -e; • V X V = (6) 
V X w = e;V • -y- (e; • V)t7= (7) 



where p is the pressure. Let me deflne the speciflc pressure, p 
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since nothing depends on z. 
Thus we can define 

w = \I^ (8) 

where i/) also obeys V^V' = and where 

Vi/i • V0 = (9) 
• = w • w = t;^ (10) 
V?/i • V-0 = w • u? = -y^ (11) 

(12) 

Let me now define a new coordinate system. I could use and V'j but I will be interested in fluid flows where the 
velocity approaches a constant value Vx = v^s^ -y^ = at large distances. I will thus instead use the functions ■)/') 
defined by 

(13) 

Vo 

V' = - (14) 

Vo 

as the new coordinates. This choice will also allow me to take the limit as the velocity vq goes to zero, where the 
potentials (/), "0 £^re undefined. Thus at large distances, (j) = x and ~ y. The spatial metric in the xy coordinates is 

ds^ = dx^ + dy^ = gijdz^dz-' (15) 

(where the Einstein summation convention has been used where a repeated index implies summation over that index, 
and where z^ = z"^ = y). Do not confuse z^ with the horizontal direction z which nothing depends on. The 
Laplacian is for a general metric function of gij{z) is 

^ ^d^^W^d, (16) 

where are the components of the matrix with is the inverse to the matrix of coefficients gij and where g is the 
determinant of the matrix with coefficients gij. For a reference regarding metrics and the coordinate independent 
equations see almost any book on General Relativity 

In two dimensions, if gij = fcjij where / is some function of the coordinates z*, then since g''-' = jg^^ and 

g — det{gij) ~ f'^det{gij) — f^g, we have = 7^^' l^otrics such as g^j and gij are said to be confomally related. 
Recalling that the change in the metric components from one coordinate system z* to a new system z^ are given by 

where the Einstein summation convention has been used, the components of the upper components of the usual flat 
space metric in this new z^ = 4>,z^ = ?/; coordinate system are 

= V(/) • V0 = 4 (18) 

g^^ ^^^j.^^^^ (19) 

g'l"!' = V0 • VV- = (20) 
le, the new metric (the inverse of this upper form metric) in these new coordinates is a conformally flat metric 



Since in the xy coordinates the metric is flat, this metric is also flat in i/j, (j) coordinates, (the curvature is not changed 
by a coordinate transformation) and the scalar curvature in this new coordinate system is zero. Using the equation 
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for the scalar curvature of a metric (and in two dimensions, the scalar curvature is the only independent component 
of the curvature) one gets 



(a| + a|)ln(4 = (22) 

(This is valid as long as ^ is not equal to zero anywhere) 
I define 

V' = dl + dl (23) 

The Laplacian 

-di^/§g'^dj^ (24) 



Via) 

is, since the metric in ip, </> coordinates is conformally flat, just 

-^VH (25) 

for any scalar function $. 

Since in x, y coordinates, the Laplacian of both the scalar functions x and y are zero, they must also be zero in cj), tjj 
coordinates ( since the Laplacian is an invariant scalar operator), and, as functions of (jijtp, we have 

V^x{(l),iP,t) =\7^y{(j),^j,t) =0 (26) 

as the equations of motion obeyed by x and y in these new coordinates. 

ip is the stream function, and the vector v is tangent to the surfaces of constant ip. v ■ Vip = v- w = 0. The bottom 
of the flow must be tangent to the flow vector (no flow can penetrate the bottom), and thus must be a surface of 
constant tp, which I will take to be -0 = 0. Similarly, if the flow is stationary, the top of the water, no matter how 
convoluted, must also lie along a streamline, since a particle of the fluid which is at the top, must flow along the top 
(the velocity of the particles must be parallel to the top surface). This means that the top of a stationary flow ( but 
not a time dependent flow) also is at a constant value of ip which I will label ipx- 

We also have 

dx<l) = = — (27) 



dy(j> = -dxip = ^d^x = d^y = (28) 



-d^x = d^y=^ (29) 



and thus 
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1 

{d^x)^ + {d^xf 



(30) 
(31) 



Solving for x and y as a function of V', <^ , which is just solving the Laplacian in terms of ijj, <l>, gives us the velocity 
at all points. 

The boundary condition along the bottom for these functions must be that the velocity along the bottom be parallel 
to the bottom. If the bottom has the functional form y = F{x) then y{(p, 0) = F{x{(f), 0)). On the top of the flow, we 
have the boundary condition that p = 0. Prom the Bernoulli equation for a stationary flow is 

+ gy + p = const (32) 

which, if the flow has constant velocity u over a constant depth bottom of height h far away from the obstacle, gives 
the equation for the top of the flow 

^v{(j),t{jTf + gy{(l),tpT) = '^vl+ gh (33) 
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Writing this in terms of (j), ip we have the upper boundary condition of 



(34) 



This is a complicated, non-linear, boundary condition. Thus while the equations of motion of x, y are simple 
(Laplacian equals zero), the physics is all contained in the boudary conditions. 

If we are given y{x) as the equation for the bottom, the solution of the above non-linear boundary value problem 
is difficult. However if , instead of specifying the lower boundary, one specifies the shape of the upper boundary 
y{(p,ijjT), one can use Bernoulli's equation in these new coordinates to determine the ip derivative of y . Since 



^1 '"v 



(35) 
(36) 



we have 



(^vy)^ + [dphiyf 



(37) 



and Bernoulli's equation is + gy = const along the top surface of the fluid where p = 0. Soving for d^y we get 

d^y{<i>,i,T) = -J-{{d'^y{<l>MY 



vl + g{y{oo, ipr) - y{(p, tpr)) 



(38) 



Any function H{ijj(l)) which is a solution of d'^^H+d'^H = can be expanded in exponentials e*'^'^. We sec immediately 
that the dependence of these modes of ip must be in terms of e^*"'' or equivalently in terms of cosh{ktp) and sinh{ktp) 
for the V dependence. Thus, since y obeys that equation, we have 



2/(^,^) = j < 



e""^ (afe cosh(/e(V; - V-t)) + sinh(fc(V' - Vt))) 



with 



= ^ j \d^y{<j>,i>T)e-"''>'d<p 



Then at the lower boundary. 



x{<l>,0) = j j 



y{k)cosh{ktpT) - d^y{k) 



smh{kipT) 



'^'t'dk 



dy{k)kcosh{kipT) + y(A;)sinh(fcVT)l e''"l'dkd(p 



where 



dy{k)= J d^y{cj>,^T)e-""^dcf> 



(39) 

(40) 
(41) 

(42) 

(43) 

(44) 
(45) 



This gives the bottom as a parametric set of functions of cp. 

In figure 1 we have an example of sub to supercritical flow over an obstacle, calculated as above. Note that the 
obstacle is a reasonable function y{x). 
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FIG. 1: Figure 1. The upper graph gives the top and bottom {y{tpT) and y{0) of a symmetric fiume flow with vo = .3m/s. 
The top of the flow was specified with y((?i,VT) = .015(e''''" -'"'' + e^'l'+-^'> Z^. Note that the bottom of the flume is a reasonable 
function of x. The lower graph gives the velocity of the fluid flow, {v{cj))) as a function of x and the phase velocity of long 
wavelength waves \/ g{y{4'^ V't) — J/(</', 0)) as a function of x. The ratio of these two velocities is the Proude number, which is 
greater than unity over the obstacle. 



A. Vo = Q limit 

The boundary condition equations are easily solved in the limit as vo — 0. The upper boundary condition becomes 
simply y = h and d^y = 0. This can be solved (in terms of the unknown lower boundary solutions 0), x{(j), 0) by 



j/(<^,v) = j afce"'^ " :" dk 



-■I 



sinh{kiljT) 



sinh{kipT) 



dk 



(46) 
(47) 



where 



"fc = ^ / y{^,o)e-""f'd<p 



(48) 



Of course, we are not given y{(j), 0) but rather y{(p, 0) = F{x{(j), 0)). However one can get rapid convergence by iteration 



xn{(j),Q) = (j) 
yi+i(<^,0) =F(a;,(0,O)) 



(49) 
(50) 
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which gives via the above equations the solution yi+i{(f), ^) and thus 

Xi+i(0,O)= d^y,+i{(j),0)d(j) (51) 



For small vq, one can get a first order correction for the surface value of y{(j),ipT) by taking 

y{<P, i,T) = h-vl -i— ^ (52) 

le, for slow flow over a bottom boundary, the stationary solution for that flow is easy to find. 

B. Formal General solution 

The general solution to the equation V^-F=0 can be written as 

F = /(<^ + i^)+5(0-iV) (53) 

If F is real, then g{(j> — iip) — {f{4> + itp))* We then have 

x{(l),ip) = xicji + ii/j) +x*{4> + itp) (54) 
y{(l),i>) = i{x{(j) + iijj) - x*{(f) + iij))) (55) 

Given the boundary conditions along the bottom, we have 

x{<i>) = ^{xo{cl),Q)-iyo{cl),Q)) (56) 
This of course still leaves the highly non-linear boundary conditions at the top to solve to find x and y everywhere. 

II. FLUCTUATIONS 

Let us assume that we have a background solution to the stationary equation, a;o(0: V'): J/o('/'j "0)1 or equiva- 
lently, ^Q{x,y),^po{x,y). We want to find the equations for small perturbations around this background flow. Let 
us also consider a solution to the full time dependent equations, (f){x,y,t),^p{x,y,t) together with their inverses, 
x{(j),tp,t),y{4>,^p,t), such that y{4>{x,y,t),t{j{x,y,t),t) = y and x{(j){x,y,t),^p{x,y,t),t) = x. Define the small de- 
viations from the background by 

(57) 

(58) 
(59) 
(60) 
(61) 

Then we have 

y = y{(j)oix,y)+6(j){x,y,t),il;Q{x,y)+Sil;{x,y,t),t) (62) 
= yo(0o(a;,2/) + S(f){x,y,t),ijjn{x.y) + 6i:{x,y,t),t) +6y{(j)o{x,y) + S(l){x,y,t),tpo{x,y) + Stp{x,y,t),t) (63) 

Keeping terms only to first order in " S" , we have 

y = yo{(l>o{x, y), il>o{x, y)) + d^yo{4>o(x, y), ipoix, y))54> + d^yo{(t>o{x, y), il>o{x, y))Sil) + 5y{(t>o{x, y),ipo{x, y)) (64) 

or 

6y{(t)o{x,y),ti)o{x,y)) = -^^5(l){x,y) - ^^^Stp{x,y) (65) 



6^ 


= (t>{x,y,t) - 


Mx,y) 




= '>P{x,y,t) - 


■^o{x,y) 


Sx 


= x{ip,(j>,t) - 


xo{(f>,ip) 


Sy 


= y{4>,'4^,t) - 


yo(0,f/') 



(where all velocity components are those in the background flow). 
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Similarly 



and 



The Bernoulli equation is 



^x = -^^^S(t){x,y) + ^^^Sip{x,y) (66) 

S<j){xo{(l),tp),yo{(l),il>)) = —{va:6x{(j),tjj) + Vy6y{(p,'tp)) (67) 
Sip{xo{(t),tp),yo{(t>,'4')) = —{-Vy5x{(f>,il)) + v^Sy{(t),'ip)) (68) 

vod,Hxi<l>, ^, t),yi<P, ^, i), + I ^d,x{c^,^,twl^d,y{^^i:,t))^ + 9y{<>^ V-, + P = const (69) 

where the first dt is defined as the derivative keeping x, y fixed, not ^, fixed. Here p is the specific pressure. 
Writing this equation perturbatively, we have 

- v^d,5x - vyd,y - + (a,^)2)2 ( V^^*'^. + + ^.Jy + = (70) 

where all of the velocities are the values of the background velocities at the location (j),^!). le, Vx{(t>,'4') = 

vox{xo{4','^),yo{(t>,'4'))- 

We can now rewrite this equation in terms of 6(p = 6(p{xo{(p, tjj),yo{(f), ip)) to get 

vodtS^ + (v^d^i^Scj) - J^V) + Vyid4,{^6cp + - + ^¥ + Sp = (71) 

RecalUng that a^ff = d^d^yo = and d^,^ = -d^^, we finally get 

vodtScf) + v^d^Scj) + 5^(iu^ + gyQ)5(f> - d^{gyo + ^v^)S^ + 5p = (72) 

The boundary conditions at the bottom are that 6x and dy must be parallel to the bottom, or VxSy — Vy6x = 
which is just 

Si;{(l),0) = (73) 

At the top, the pressure at the surface must be 0. However the surface is no longer simply ^ = i/jt because of the 
time dependence of the equations. Let us assume that the surface is defined by 

V' = ^'(0,t)+V'T (74) 

Since a particle of the fluid which starts on the surface, remains on the surface, we can define the fiuid coordinates 
r],(- Then the velocity of the fluid is 

= j^cl>{C,v,t) (75) 
v'^ = j^HC,V,t) (76) 



Along the surface, we therefore have 



But, 



= dt'i' + v'>'d^^ (77) 



v'^ = j^<t>{x{rj,C,t),y{v,C,t),t) (78) 
= Vxdx(l> + Vydycj) + dtcj) (79) 

= —+dt4>{x,y,t) (80) 
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Thus, assuming that ^ is also small (the same order as the other "5" terms), we have 



(81) 



vo 

On the surface, we have the Bernoulli equation, which to first order is 

^v^{<P, V-T + *) + 9yo{<P, V-T + *) - l^'^i'P, + *) + 9yo{<P, V-t + *) + dtS(j) 
+ v^d^5(j) - d^{^v'^ + gy)5(f) - d^i^v^ + gy)5tp + p - pq = 

But along the surface = -(/t, the background ^v'^ + gy is constant, so the (f) derivative is 0. We have 

{dt + v''d^)5cp + d^i^v^ + gy){^ - ^^) = 
Dividing by G = d^{^v'^ + gy) and taking the derivative dt + we get 



Vo 



1 ~ w2 

G Vo 



—d^Si> = 

Vo 



as the equation of motion for the surface wave. 6(p and 6ip are related by the boundary condition 
bottom. 

Since both 5(f) and Sip obey V'^Sip = V'^Scf) = 0, we have 



Furthermore, since 



so 



V^5V = V^5?!> = 



dyS(f) = —OxSip 



d^Scj) = d^xodxSip + d4,yody54>(t> 

= d^yodydip - d^xoi-dxSip) = d^Sip 
d^5(f) = -d^6il> 



(82) 

(83) 
(84) 

(85) 

(86) 
along the 

(87) 



(89) 



(90) 
(91) 
(92) 



For irrotational time- independent flow, the acceleration of a parcel of fluid is v ■ Vv = y{^v ) and the orthogonal 



component of this, the centripetal acceleration is 

1 



(93) 



Also gd^y = g^^^ « gvo/v soGv/vq is the effective gravitational field orthogonal to the flow lines (including the 

centripital acceleration) . 

However it is important to note that it is the effective force of gravity only at the surface of the fluid, not at the 
obstacle to the flow along the bottom, that is important for the equations of motion. 



III. SHALLOW WATER WAVES 

Since </>, tp are real functions, the solutions can be written as 



d(l){(j), xP) = + iV) + iZ{<j) + ii'))* 
5,p{<t>, iP) = i{Z{ct> + i^) - {Z{<p + li,))*) 



(94) 
(95) 
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for some function Z. These functions clearly satisfy the Laplacian equation for, and furthermore also satisfy the 
differential relations on the derivatives oi x,y with respect to c^^-psi This gives 

= 5V(<^,0) = i(Z(</.)-Z*(,^)) (96) 

le, Z is a real function of a real arguments, which gives 

b<^{^, V) = W + ^V-) + - iV')) « 2Z(0) + Z"(<^)i? (97) 

= 2V'^'(<^) (98) 

or, to first order in t\)T 

= i)Td^b<^ (99) 

The equation for the waves then becomes 

iPt + v^d^)^ {Bt +v^d^)6^- v^i)Tdl5(t> = (100) 

We note that this is not a Hermitian operator acting on 5(f). Recall that a Hermitian operator is one such that 

j 54H(l)d(l)dt = j {'HS^)5(l)d(l)dt (101) 

if we assume that all of the boundary terms in the integration by parts are zero. We can rewrite the equation for S(f) 
by dividing by v"^ as 

{dt + d4,v^)^idt + v^d^)5<i) - ^Tdpcl) = (102) 

This is a symmetric equation, derivable from an action, 

1 



/ 



d<pdt (103) 



This action has the global symmetry (50 — >■ e^^Scf) and thus has the usual Noether current associated with this 
symmetry. In particular it has the conserved norm 

<,50,50'>=^ J i^S(l>*-^{dt + v^d^)^-^-d^'-^{dt+v^d^)^-^^d<l> (104) 

IV. DEEP WATER WAVES 

For deep water waves, we can assume that cither Z{(j> + iipr) » Z{ip — i(j)T) or Z{4) + zi/'t) << Z^tj} — i<j)T)- (ie, 
we assume that as analytic functions, Z goes to zero either in the upper or lower half plane.) 

Let us also assume it is the first case, and let us define Z{<t)) = Z{(j) + iipT), and that dtS(p = iui6(j) We then have 

{iuj + v''d4,)^{iiJ + vl)]Z - i-iyd^Z = (105) 

If we assume that K = i{d^ln{Z)) is large and negative, such that Z varies faster than v'^ or G, we have approximately 

G 

or 



+ Kv^ = (106) 



oj = -v'^K ± Vv'^GK (107) 
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V. GENERAL LINEARIZED WAVES 

The equation in general is 

{Bt + v^d^)-^{dt + v'^d4,)5(t> - v^d4,6^ = (108) 

Fourier transforming with respect to (p and ip, and using the fact that ^■^ = at = 0, the functions 6(p, 6tp then can 
be written as 

5(j){(j),i},t) = j A{k,t)e'^'t' cos'h{k'4})dk (109) 

di}{(j),'ip,t) ( A{k,t)e^'"^smh{ki})dk (110) 



since again they obey the Laplacian equal to zero in these variables. 
Defining B{k, t) = A{k, t)cosh{ktp) this can be written as 



5(t>{^,tl>T,t) = J B{k,t)e''">'dk (111) 

5^{^,tpT,t) = i J B{k,t)e''"l'^ tainh{ki))dk = ita.nh{-ii;Td^) j B{k,t)e''"l'dk = iteinh{-i^Td4,)6(() (112) 
Thus the equation of the surface waves can be written as 

= {dt + d^v^)-^{dt+v^d4,)S(t)-idkt&nh{-ii}Td^)5(l> (113) 

le, we get the usual tanh dispersion relation for the transition from shallow to deep water waves. 
This equation is symmetric and real, and thus if S(j) is a solution, so is 6(p*. Again this gives a conserved norm 
between two solutions to the equations of motion 5(j) and Sc/)' of 

< S(t>, 5(1)' >= J ]^(j>*{dt + v^d4,)6<j,' - 5ci>'{dt + v''d^)5(l)*'\ dcj) (114) 

We note that this equation depends only the conditions at the surface of the flow. It is defined entirely in terms of 
the factors [v"^ and G = d^{gy + ^v"^) defined at V' = i'T, and is independent of the obstacles, or the fiow throughout 
the rest of the stream except insofar as they affect the flow at the surface. This might well change if either vorticity 
or viscosity were introduced into the equations. 

This norm is crucial to the analysis of the wave equation. It is conserved (in the absence of viscosity), and in 
the use of such waves as models for black holes, it is this norm which determines the Bugoliubov coefficients (or the 
amplification factor) for waves in the vicinity of a horizon (blocking flow in the hydrodynamics sense) and determines 
the quantum noise (Hawking radiation) emitted by such a horizon analog. The quantum norm used in the quantization 
procedure is 

{S<p,6<l>)Q='-<dcP,5cP> (115) 
If we deflne a new coordinate 4> = J ^dcj), the norm becomes 

< 5cf>, 5(t>' >= j (^ [6(l>*{dt - d^)5cl>' - 5cj>*{dt - d^)6cf^ # (116) 
If the surface of the fiow is shallow (^ << 1) then ^ = Vx ~ v and 4>^ 

To relate this to the measured quantity, the vertical displacement at the surface of the waves, we must relate 5(j) to 
5y at the surface of the fluid. We have 

*(i, (j)) = il){t, X, yrit, x)) - Vt = Stp{t, x{(i), Vt), 2/(0, tpr) + Vxdyr (117) 

or 

Syr = — (* - ^V) = -T^dt + v^d^)5(l> (118) 
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Now, Gvx ~ ff^ ~ 5 (ignoring the centrifugal contribution to the effective gravity), so the norm becomes 



< 5yt,5yT > 



{dt + dr)-^Sy*rSyT - {dt + dr)-^SyTdy*T 



{{dt + d.)-'V^Sy^)V^SyT - {{dt + d.)-'V^SyT)V^Sy^] d4> 



(119) 
(120) 



and d^^ « ^ 

If we assume that the incoming wave is at a set frequency lj and take the fourier transform with respect to t, x of 
^ v{ip)yT{t, 4>) this becomes 



< 5y,5y >-- 



l(V^yT)(fc)P 

{uj + k) 



dk 



We can also look at the norm current. 



V 



(121) 



(122) 



/ 

J6 



dA^{dt + v^d;)6(t> - d^{^{dt + v^d^)S^*) + [(-i5^tanh(-i^T5^)^<^*)^<^ - (i9^tanh(#T)a0)^<^*)^<^] # 



The integrand is a complete derivatives. Although this is not obvious for the terms with the tanh in them, we can use 

2n-l 



{df^sm - 5rdf6<i> = J2 (-i)'-aj5<A-a|"-i-'-^0) 



(123) 



r=0 



and the fact that id^ ta,nh{itpTd^) can be expanded in a power series in 9| to show that they also a complete derivative.. 

Thus the integrand can be written in terms of a complete derivative of with respect to and we can regard the 
term that is being taken the derivative of as a spatial norm current J'^ so that if J* is the temporal part of the norm 
current, we have 9/ J* + 9^ J"^ = 0. 

If we are in a regime where d(j) = ^e"'"*"*^*^, (ie, a regime where the velocity v and G are both constants), then we 
have 



Gv^ Gv^ 



(124) 



where Vp and Vg arc the phase and groiip velocity of the wave. In a situation in which one has a wave train with 
some definite frequency and wave number entering a region, then the sum of all the norm currents for each k at the 
boundary of the region must be zero. 



VI. BLOCKING FLOW 



Let us return to the static situation. Define U 



we have the equation 



d,p--^U + itaTih.{i'il;Td,j,)U = 
G 

As above, there is a solution if we assume that the derivatives are small, which gives 

const 



U 



For rapid variations, we have 



i7 = const^e^/^''* 



(125) 



(126) 



(127) 
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with the transition from one to the other occuring roughly when the logarithmic derivatives of the two solutions are 
equal 



^'''^ 

(129) 

Defining the Froude number by = (the square of the velocity of the fluid over the velocity of the long 

wavelengths in the fluid in the WKB approximation), we have 



- ■ '4(ln(F)')2 + ] (130) 



F2 - 1 V ^ ' ^ ' \F^^j 

Note that for a non-trivial rate of change of of the bottom, the turning point occurs well before the horizon. 

The ' denotes derivative with respect to not x. We can rewrite this approximately (assuming that ^ « 1 and 

that 2^"(^) <^ and tpr ~ vd where d is the depth of the water at postion x. as 



dF^ _ {F^ - 1) 
dx ~ F'^d 



(131) 



Note that this transition occurs before Gipr = or Froude number equals 1. The wave on the slope piles up and 
its frequency makes the transition to deep water wave before we hit the effective horizon. 
The long wavelength equation, 

^ {dt + v^d^)^ {dt +v^d^)6ct>- ^TdlScl> = (132) 

is not that of a two dimension metric, which is always conformally flat, but can be written as a the wave equation for 
a three dimensional metric where all derivatives are equal to zero in the third ^ dimension for the variable 6(j). The 
metric is 

ds^ = a((l - -^)dt^ + 2-^dtd<l> - ^7^#') - -^^de (133) 

where a is an arbitrary function of (j), a two dimensional conformal factor which does not affect the two dimensional 
wave equation This metric has surface gravity 

^, .\ 1..2« Z.2 



(The surface gravity is the acceleration in the metric as seen from far away, for a static time independent metric in a 
coordinate system which is regular across the horizon, it can be deflned by k = r|j at the horizon, where Vj^. is the 
Christofell symbol for the metric. Then r*j = — \g*''^{d^gtt) at the horizon.) 

VII. CONVERSION TO 5y 

Of course is not what is actually measured in an experiment. That is the fluctuation l^y{x) which is the difference 
in height between the stationary flow, and the height with the wave present. We can relate this to (5^' and 

ys{x, t) = yo{x) + Ay{x, t) (135) 

where yo is the surface for the background. 

6y = \5(l> + %5i! (136) 

Since 5il} = t&uh{'^ Hd^)5(j), we have 

v^6y = [vy + Hd^))]6(t) (137) 
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Inverting this for deep water waves, 



while for shallow water waves 



Vy+Vx 



Sy (138) 



exp^ ^'^^—Sydcj) (139) 



The integrand in the exponent is non-zero only in the region where the background flow is dimpled, and, since ^ is 
in general very small, the exponential can be neglected in most situations. 

In the intermediate region, where the wave changes from shallow to deep water wave, there is no easy solution to 
these equations, but they can be integrated numerically. 



VIII. WAVES IN STATIONARY WATER OVER UNEVEN BOTTOM 

In the limit as vq goes to zero, so does v with the ratio being a finite function, y obeys the equation + d^y = 
with the boundary conditions along the bottom that y ~ Y{x). with Y the given function of x of the bottom, and 
along the top, y = H, a, constant. If we assume that we know Y{(p) (instead of Y{x)) along the bottom, this can be 
solved by 

+ (140) 



where 
and 



a{k) = ^ I Y{4>)e''"^ (141) 



One gets rapid convergence if one starts by taking x = (j), substituting into Y{x{(j))) to find Y{(j)), finding the new 
x{(j)) and substituting in again. 

Then ^ at the surface is zero, while 

The equation for small perturbations becomes 

- id^ tanh(iVTa^)<5</) = (144) 

where 

-^G = —gd^y = g— = g-^ (145) 

.If the depth is constant, the backgound = y and (f) = x giving the usual equation, which allows us to write 

d'i5(t) + igdxtsjiH:4jTd4,)5(f) (146) 

For deep water waves, where the tanh is unity, this equation is exactly the same as the deep water equation for 
constant depth. The fact that the bottom varies makes no difference to the propagation of the waves, as one would 
expect. 
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For shallow water waves, where the tanh can be approximated as the linear function in its argument, the equation 
becomes 

5^50 = xPrd^d^Si^ = g^j,'!^dl5c^ (147) 

This allows us to determine the wave propagation over an arbitrarily defined bottom. Note that in the stationary 
limit, the background flow is certainly irrotational, implying that the assumptions made here should certainly be valid 
(of course neglecting the viscosity of the fluid). 
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